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ABSTRACT  7 

^The  effect  of  dust  on  aerodynamic  drag  of  spheres  and  cylinders 
was  calculated  by  using  a mathematical  model  developed  for  this  purpose. 
The  results  were  compared  to  experiments  previously  done  by  other  workers. 
The  calculated  and  experimental  results  agree  favourably,  showing  that 
the  mathematical  model  is  satisfactory.  Impaction  efficiencies  and 
drag  coefficients  due  to  dust  alone  were  then  obtained  using  the  model 
for  a wide  range  of  the  Inertia  parameter  and  the  results  are  presented 
graphically.  The  model  can  also  be  used  for  calculating  velocity  dis- 
tributions and  points  of  Impact  for  a stream  of  airborne  particles  flowing 
over  a sphere  or  cylinder. 
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frontal  area  of  target,  cm2 

far  upstream  cross  sectional  area  of  the  envelope  of  particles 
which  eventually  hit  the  sphere  or  cylinder,  cm2 

tube  radius,  cm 

dimensionless  drag  coefficient  for  spheres  or  cylinders  In  air 

drag  coefficient  of  cylinder  or  sphere  due  to  stream  of  particles 
alone 


total  effective  drag  coefficient  for  spheres  or  cylinders  due 
to  particles  and  air  together 

drag  coefficient  for  spheres  in  tubes  adjusted  to  a sphere  In 
free  space  by  the  Landenburg  correction 

drag  coefficient  for  particles  alone  adjusted  from  v to  U 

particle  diameter,  cm 

Impaction  efficiency  of  particles  on  a target 

force  exerted  on  a target  by  the  air  and  total  stream  of  particles 
hitting  It,  dynes 

force  exerted  on  a target  by  the  air  alone,  dynes 

force  exerted  on  the  sphere  or  cylinder  by  the  total  stream  of 
particles  which  hit  It,  dynes 

force  exerted  on  a target  by  the  particles  alone  travelling  In 
air,  dynes 

cylinder  or  sphere  radius,  cm 

total  mass  of  particles  which  had  Its  momentum  changed  In  unit 
time,  g s"1 

drag  pressure  on  sphere  or  cylinder  due  to  stream  of  particles 
alone,  dynes  cm"2 

distance  from  any  point  to  the  origin,  cm 
time,  seconds 

free-stream  velocity,  cm  s"1 
local  fluid  velocity,  cm  s_1 
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radial  component  of  fluid  velocity,  cm  s"1  ■[ 
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local  particle  velocity,  cm  s- 1 | 
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rms 

root-mean-square  free  stream  particle  speed  upstream  near  target  1 

position,  cm  s_1  f 
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vx 

component  of  particle  velocity  parallel  to  free  field  flow  | 

direction  Immediately  before  Impact,  cm  s"1 
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component  of  particle  velocity  parallel  to  flow  direction  1 

Immediately  after  Impact  and  reflection,  cm  s"1  | 
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far  upstream  particle  velocity,  cm  s"1  ii 
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vy 

component  of  particle  velocity  perpendicular  to  free  field  i 

flow  direction  Immediately  before  Impact,  cm  s"1  1 
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Avx 

change  In  the  component  of  particle  velocity  parallel  to  the  1 
free  stream  flow  direction  caused  by  reflection  from  the  sphere  I 
or  cylinder,  cm  s'1  l 
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X 

co-ordinate  (origin  at  centre  of  sphere  or  cylinder)  of  if 

particle  position  parallel  to  free  stream  flow  direction,  cm  1 
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transverse  co-ordinate  of  particle  position,  cm  f 
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off-axis  distance  of  a particle  at  point  of  Impact  with  the  .) 
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far  upstream  transverse  co-ordinate  of  the  envelope  of  particles  j 
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angle  of  Incidence  of  a particle  against  the  sphere  or  cylinder, 

radians  ! 
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angle  between  particle  direction  and  free  field  flow  direction 
just  before  the  particle  hits  the  sphere  or  cylinder,  radians 
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total  angle  change  of  a particle  from  Its  far  upstream  direction  to 

Its  direction  after  reflection  from  the  cylinder  or  sphere,  radians  ' 

1 

0 

polar  angle  between  x axis  and  radius  vector  to  particle  position, 
radians  ) 

U 

absolute  velocity  of  fluid,  poise 
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fluid  density,  g cm"3 
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total  mass  of  particles  and  air  per  unit  volume  of  space,  g cm-3 
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uniform  particle  density  per  unit  volume  of  air  far  upstream, 
g cm"3 
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particle  density,  g cm-3 
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The  following  were  dimensionless: 


f/-,  change  In  component  of  particle  velocity  due  to  reflection  from 

'y'  sphere  or  cylinder  as  a function  of  off-axis  distance  far  upstream 

K particle  Inertia  parameter 

Re  spherical  particle  Reynolds  number  In  flow  Influenced  by  presence 

of  sphere  or  cylinder 

ReQ  spherical  particle  Reynolds  number  In  free  stream 
r , distance  from  any  point  t:o  the  origin 

T time 

u local  fluid  velocity  jj  , ux  and  uy  are  x and  y components 

v local  particle  velocity  jj  , vx  and  vy  are  x and  y components 

vx  gy  , parallel  component  of  particle  velocity 

Vy  , transverse  component  of  particle  velocity 

x parallel  co-ordinate  y 

y transverse  co-ordinate  jj 

yd  far  upstream  transverse  co-ordinate  of  the  envelope  of  particles 

which  eventually  hit  the  cylinder 

y,  co-ordinate  of  particle  used  In  discretization  for  Integration; 

y.|  Increases  with  1 from  y^  ■ o to  yn  * yd 

4>  dimensionless  group  independent  of  drop  size  formed  by  combining 

ReQ  and  K 
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1.  INTRODUCTION 

Extensive  studies  have  been  made  of  the  drag  exerted  on  spheres 
and  transverse  cylinders  by  moving  air  (Schllctlng,  1960)  and  attempts 
have  been  made  to  determine  the  effect  of  introducing  dust  Into  the  air 
stream  (Gillespie  and  Gunter,  1957;  1959).  Also,  the  trajectories  of 
water  drops  moving  In  the  neighborhood  of  a circular  cylinder  placed  In  a 
uniform  stream  of  air  have  been  calculated  by  numerical  solution  of  the 
equations  of  motion  (Glauert,  1940)  and  with  a mechanical  analog  (Brun 
et  al.,  1953).  The  latter  method  has  been  used  in  conjunction  with  flight 
Instruments  used  to  study  droplet  size  and  distribution  In  Icing  clouds 
(Brun  and  Mergler,  1953).  Measurements  of  aerodynamic  drag  on  circular 
cylinders  due  to  the  blast  wave  from  large  TNT  bursts  have  been  done  at 
DRES  (Mel  1 sen,  1974).  The  results  of  these  tests  have  been  used  to  prpvlde 
blast  loading  Input  for  structural  analysis  of  lattice  type  masts  also 
tested  In  the  trials.  On  some  of  the  trials  there  have  been  visible 
quantities  of  dust  associated  with  the  blast  wave.  In  general  the  measured 
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strains  in  the  masts  have  agreed  favorably  in  an  Individual  field  trial 
with  the  ones  predicted  from  theoretical  analysis  based  upon  the  blast 
loading  Input  from  the  drag  data  averaged  over  several  blast  trials 
(Laidlaw,  1977).  The  Influence  of  dust  Is  difficult  to  determine  and 
generally  considered  negligible. 

Another  area  of  Interest  In  connection  with  moving  particles 
In  an  air  stream  Is  the  Impingement  of  liquid  drops  on  protective  clothing 
and  military  equipment.  The  Impaction  efficiency  and  force,  and  the 
droplet  distribution  are  all  of  Interest.  The  study  of  these  effects  on 
transverse  cylinders  and  spheres  Is  useful  In  this  application. 

The  purpose  of  this  report  Is  to  describe  numerical  solutions 
to  the  equations  of  motion  of  particles  In  a moving  air  stream,  with 
possible  applications  to  both  the  effect  on  aerodynamic  drag  from  the 
point  of  view  of  blast  loading  of  structures  and  the  Interaction  of 
droplets  with  spheres  and  cylinders  from  a chemical  defence  point  of  view. 
The  mathematical  techniques  along  with  their  associated  computer  pro- 
grams are  described  herein.  Also,  the  results  are  compared  with  experi- 
ments done  at  ORES  (Gillespie  and  Gunter,  1957). 


2.  DEFINITION  OF  THE  PROBLEM 

A spherical  particle  flowing  at  free  stream  velocity  far  up- 
stream of  a spherical  or  transverse  cylindrical  solid  target  will  not 
necessarily  follow  a streamline  In  the  vicinity  of  the  target  where  the 
radial  and  axial  velocity  components  of  the  fluid  may  be  changing  markedly. 

A particle  that  just  Impinges  on  the  outer  circumference  of  the  target  was 
considered.  All  particles  of  the  same  diameter  within  an  envelope  generated 
by  this  limiting  particle  trajectory  would  hit  the  target,  assuming  that  no 
particles  collide  with  other  particles  reflected  from  the  target. 

All  particles  In  the  steady  flow  field  far  upstream  were 
assumed  to  be  travelling  at  the  velocity  of  the  fluid  but  the  velocities 
of  the  particles  will  differ  from  the  fluid  velocity  as  the  particles 
approach  the  target.  This  Is  due  to  Inertial  effects  of  the  particles  In 
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the  flow  field  near  the  object.  The  velocity  of  a particle  immediately 
before  contact,  and  also  the  point  of  contact  with  the  target,  depend 
upon  the  position  of  the  particle  within  the  limiting  envelope  far  up- 
stream from  which  It  travels  to  the  target.  These  particles  all  have 
their  momentum  changed  when  they  collide  with  the  target  and  in  the  case 
of  solid  particles  are  deflected  In  a new  direction.  The  total  change  In 
momentum  Is  balanced  by  the  drag  force  exerted  by  the  particles  on  the 
target. 

The  basic  problem  was,  therefore,  to  determine  the  envelope  of 
particles  hitting  the  target  and  the  change  In  particle  motion  due  to  the 
presence  of  the  target  In  the  flow.  Then  the  drag  coefficient  due  to 
particles  alone  could  be  obtained. 


3.  EQUATIONS  OF  MOTION 


The  motion  of  an  Individual  spherical  particle  has  been  shown 
(Batchelor,  1967)  to  be  determined  by  the  following  ordinary  differential 
equations: 

tiy  m ye(uv  - \)  (Eq.  ,) 

dT  24K 


where 


^x  . C0Re<G»  • 5x> 
dT  24K 

Re  = ReQ[(uy  - vy)2  + (ux  - v^2]*4 

j 2i  i 

K » - — particle  inertia  parameter 
18yL 

Re0  « free  stream  Reynolds  number 


(Eq.  2) 

(Eq.  3) 
(Eq.  4) 

(Eq.  5) 


The  symbols  are  defined  In  the  notation  section  near  the  front 
of  this  report  and  the  basic  geometry  of  the  flow  system  Is  Illustrated 
in  Fig.  1. 
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Several  assumptions  were  Inherent  In  the  development  of  Eq.  1 
and  Eq.  2 for  calculating  Impaction  efficiency  and  force  due  to  a stream 
of  particles.  Including. 

(a)  uniform  particle  distribution, 

(b)  no  gravitational  or  electrostatic  forces  of  consequence, 

(c)  monodlsperse  spherical  particles  with  diameter  very  snail 
In  relation  to  the  target  sphere  or  cylinder  diameter,  and 

(d)  free  stream  flow  that  was  steady.  Incompressible  and 
Irrotatlonal . 

The  drag  coefficient  Is  a function  of  Re>  jlds  number  and  was 
available  In  the  form  of  definitive  empirical  equations  (Davies,  1945). 
These  equations  are  stated  as  follows: 


CDRe* 


- 2.3363  x 10"4(CDRe2)2  + 2.0154  x 10"6(CDRe2)3  - 6.9105  x 


10“9(C0Re2)4 


(Eq.  6) 


for  Re  < 4 or  CpRe2  < 140 


logjoRe  ■ -1.29536  + 9.86  x 10” 1 (log10CDRe2)  - 4.6677  x 10“2(log10CDRe2)3  + 
1.1235  x 10’3  (logioCpRe2)3  (Eq.  7) 

for  3 < Re  < 104  or  100  < CQRe2  < 4.5  x 107 

4.  AIR  FLOW  FIELD 

The  assumption  was  made  that  the  air  flow  In  front  of  the 
cylinder  or  sphere  Is  given  by  classical  hydrodynamics  theory. 

a)  Cylinder 

The  equations  of  fluid  velocity  were  derived  from  the  stream 
function  for  Ideal  flow  around  a cylinder  (Mllne-Thomson,  1960).  These 
were  normalized  to  the  free  stream  velocity  and  cylinder  diameter  and 
written  as  follows: 


ux  ’ 1 


x2  - y2 


(Eq.  8) 


u - 1 - 


(Eq.  9) 
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b)  Sphere 

For  purposes  of  comparison  of  the  results  of  the  theoretical 
work  reported  herein  to  existing  experimental  data  the  flow  field  about 
a sphere  on  the  axis  of  a circular  tube  was  applied.  The  flow  about  a 
sphere  In  the  absence  of  the  tube  was  then  readily  available  from  this 
as  a limiting  case. 


The  flow  velocity  was  derived  (Mellsen  et  VI.,  1966)  from  the 
vector  potential  (Smythe,  1964).  The  velocity  components  are  given  In 
terms  of  a series  solution.  For  a sphere  diameter  of  0.8  cm  and  tube 
diameter  of  2.54  cm  used  In  the  experiments  (Gillespie  and  Gunter.  1957) 
sufficient  accuracy  (Mellsen  et  al.»  1966)  was  obtained  by  retaining  three 
terms  of  the  series. 


The  equations  were 


ur 

— ■ - cos  0 
U 


1 - 


2C 


o C 3 


2A_ 


u. 


■ sin  0 


3 r a 
C.  r 3 2A. 


1 +~2.  £ 3 


(Eq.  10) 

(Eq.  11) 


The  constants  CQ,  C and  AQ  which  are  dependent  upon  the  sphere 
and  tube  diameters  (Smythe,  1964)  are  given  as  follows: 


CQ  - 1.5401075 

C “1.0  (normalized  sphere  radius) 
Irh  I(2)C_ 

A0  - -C  £ ^ ^ , 1(2)  “ 7.5098907 

a ■ C (normalized  tube  radius) 
u0  and  ur  are  shown  In  Fig.  2. 


Converting  to  rectangular  co-ordinates  gives 
u uB 

u e - — — cos  © + — sin  0 (Eq.  12) 

* U U 
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u ufi 

Uu  - — sin  e + — cos  q (Eq.  13) 

y U U 


The  equation  for  flow  around  a sphere  In  the  absence  of  the 
tube  Is  given  by  letting  the  value  of  "a"  approach  Infinity. 

For  a -*■  » , CQ  -*■  1.5  (Smythe,  1964) 


r 

Eq.  10  becomes  — 
U 

ue 

Eq.  11  becomes  — 
U 


(Eq.  14) 
(Eq.  15) 


Eqs.  14  and  15  are  well  known  In  Ideal  flow  theory  (Batchelor, 
1967)  and  therefore  serve  as  a partial  check  on  the  correctness  of  Eqs. 
10  and  11. 


5.  VELOCITY  CHANGE  OF  AN  INDIVIDUAL  PARTICLE 

To  calculate  the  drag  force  due  to  airborne  particles,  the 
velocity  change  of  an  Individual  particle  due  to  the  presence  of  the 
cylinder  or  sphere  In  the  flow  was  first  considered.  The  velocity  change 
determined  applied  to  both  the  cylindrical  and  spherical  targets  because 
the  same  geometry  could  be  used  for  both  (Fig.  3). 

The  motion  considered  concerned  a spherical  particle  which  starts 
far  upstream  within  the  envelope  of  particles  which  travel  to  the  target, 
eventually  hit  it  and  are  reflected  by  It  In  a new  direction.  The  total 
change  in  direction  of  a particle  which  comes  from  a far  upstream  off- 
axis  position  y,  and  hits  the  target  at  a point  defined  by  y'  (Fig.  3), 

Is  represented  by  y.  For  a spherical  particle  the  angle  of  Incidence, 
a,  Is  equal  to  the  angle  of  reflection.  Then 


where 


+ a 


(Eq.  16) 


0 » sin"1 


tan-1 


yl  ] 

Vl2  - (y1)2/ 


(Eq.  17) 
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and 

where 


a + 3 


3 - tan"1  r^- 


\ 


(Eq.  18) 
(Eq.  19) 


and  vx  and  vy  are  the  components  of  particle  velocity  just  before  Impact 
with  the  target.  Substituting  the  value  of  a obtained  from  Eq.  18  Into 
Eq.  16  gives 

Y * it  - 2©  - 6 (Eq.  20) 

Substituting  Eqs.  17  and  19  Into  Eq.  20  gives  an  equation  of  y suitable 
for  computation.  Thus: 


y ■ it  - 2 tan"1 
or  In  non-dimensional  terms 


Y ■ ir  - 2 tan"1 


_ 512 


- tan"1 


vr 


- tan-1 


V 

vx 


V 

v 

\ x / 


(Eq.  21) 


(Eq.  21a) 


V ' ■ V COS  Y 


The  component  of  velocity  In  the  direction  of  free  field  flow  Immediately 
after  particle  reflection  Is 

(Eq.  22) 

(Eq.  23) 

* v cos  y (Eq.  22a) 

(Eq.  23a) 


v -Vvx2  + vy2 


where 

or  In  non-dimensional  terms  v 
and 


V ■ V V,2  + Vy2 


The  total  change  of  the  component  of  velocity  In  the  direction 
of  free  stream  flow  due  to  reflection  from  the  target  Is  then  given  by 
the  following  relationship: 

avx  ■ vx  - v cos  y (Eq.  24) 


6.  THE  DRAG  FORCE  DUE  TO  PARTICLES  ALONE 

The  force  Fp  of  a stream  of  particles  hitting  a cylinder  or 
sphere  Is  related  to  the  change  In  momentum  which  In  general  is  given 
by  the  following  relationship: 
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Fp  * m'AV)<  (Eq.  25) 

where  m'  Is  the  total  mass  of  particles  having  Its  momentum  changed 

per  unit  time  and  Avx  Is  Its  change  of  velocity  In  the  free  stream  flow 

direction.  , 

The  quantity  m1  can  be  calculated  using  the  following  equation: 
m*  “ PpAdvX0  (Eq.  26) 

where  pp  Is  the  uniform  particle  density  per  unit  volume  of  air  far  up- 

stream, Ad  Is  the  far  upstream  cross  sectional  area  of  the  envelope  of 
particles  which  eventually  hits  the  sphere  or  cylinder,  and  vXQ  Is  the 
far  upstream  particle  velocity  which  Is  assumed  equal  to  the  fluid  velo- 
city. Combining  Eqs.  25  and  26  gives 

F ■ “pVxo4vx  <E*  Z7> 

avx  was  not  constant  but  was  shown  In  the  previous  section  to  be  a function 
of  the  off-axis,  upstream  starting  location  of  an  Individual  particle.  This 
can  be  accounted  for  by  Integrating  the  velocity  change  over  the  range  of 
starting  locations  across  the  particle  envelope.  Then  the  equation  of 
force  becomes  Aj 

F - ppVX0J  AvxdA  (Eq.  28) 

o 

where  dA  Is  an  element  of  cross  sectional  area  at  location  y. 


Assuming  that  the  particles  do  not  Interact  with  each  other, 
the  equation  of  drag  force  due  to  particles  alone  Is  then  obtained  In 
usable  form  by  substituting  Eq.  24  Into  Eq.  28  which  gives  the  following 
equation:  . 


Vxo 


(v  - v cos  y)d A 

A 


(Eq.  29) 


7.  THE  DRAG  COEFFICIENT 
(a)  Cylinder 

For  a transverse  cylinder  of  unit  length,  Eq.  29  becomes 
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D ■ -A 
Ps  irL2 


(Eq.  37) 


Substituting  Eq.  36  Into  Eq.  37  and  the  resulting  equation  Into  Eq.  32 
gives  the  drag  coefficient  for  a sphere  as  follows: 


4vyfl  I (vtf  - v'  cos  Y>y  dy 


(Eq.  38) 


8.  INTEGRATION  METHOD  FOR  THE  DRAG  COEFFICIENT  EQUATIONS 
a)  Cylinder 

Let  f(y)  ■ vx  - v cos  y (Eq.  39) 

Then  Eq.  33  can  be  written 

h 


* 25xo f’O 


f(y)  dy 


(Eq.  40) 


M 

? 


f(y)  dy  was  Integrated  numerically  by  dividing  the  upstream 


cross  sectional  area  Into  strips  of  equal  width  ay,  and  calculating  the 
contribution  from  each  strip.  These  contributions  were  then  summed  as 
follows: 

yd  £/  “ \ x jt/r.  \ 1*n-l 


yd  f(yJ  + f(y„) 

f f(y)  dy  - ^ f^JAy  (Eq.  41) 


1-2 


where  y^  Increases  with  1 from  y'  ■ o to  yn  ■ yd< 
b)  Sphere 


Using  Eq.  39  In  Eq.  38  gives  the  following  equation  for  the 
drag  coefficient  of  a sphere 


Cr,  " 4v 
Ds  xo 


J f(y)  y dy 


(Eq.  42) 
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f(y)  y dy  was 


Integrated  numerically  by  dividing  the  upstream 


cross  sectional  area  Into  concentric  annuli  of  equal  width  ay  and  cal- 
culating the  contribution  from  each  ring.  These  contributions  were  then 
summed  as  follows: 

f f(y)  dy  - £ 1 f ) y1  ay  (Eq.  43) 

o 1-1 

where  y^  Increases  with  1 from  y^  « ay  to  yn  ■ yd 


9.  SOLUTION  OF  THE  DRAG  COEFFICIENT  EQUATIONS 

To  solve  Eq.  34  or  Eq.  38  for  the  drag  coefficient  It  was  neces- 
sary to  know  the  values  of  vx,  v*  and  y In  the  range  0 s y s yd. 

The  first  step  was  to  find  the  value  of  yd  which  was  the  upper 
limit  of  the  Integral.  This  was  done  by  an  Iterative  procedure  called 
the  half  Interval  method  (Carnahan  et  al.,  1969).  The  value  of  y for  the 
critical  particle  was  estimated  far  upstream  and  the  path  followed  to  the 
target.  The  difference  between  the  ordinate  of  the  target  path  and  the 
ordinate  of  the  point  on  the  actual  path  parallel  to  the  tangent  path  was 
the  miss  criterion  used.  The  direction  of  the  tangent  path  was  not  known 
a priori  but  was  assumed  parallel  to  the  actual  path.  The  half  Interval 
method  previously  mentioned  was  applied  to  determine  a better  Initial 
estimate.  Then  the  path  was  followed  to  the  target  again  for  another 
calculation  of  miss  distance.  This  process  was  repeated  several  times 
until  sufficient  accuracy  was  achieved.  The  plane  of  Initial  position 
of  particles  which  was  perpendicular  to  the  flow  direction  was  located 
far  enough  from  the  target  so  that  free  stream  conditions  prevailed.  A 
distance  of  five  target  radii  upstream  of  the  target  centre  was  considered 
adequate  (Batchelor,  1967). 

The  path  of  an  individual  particle  was  determined  step  by  step 
by  applying  a fourth  order  Runge-Kutta  method  (Carnahan  et  al.,  1969)  to 
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the  equations  of  motion  (Eqs.  1 and  2).  The  values  of  Re  and  K In  these 
equations  were  easily  determined  for  each  new  step  by  direct  substitution 
of  previously  determined  values  Into  Eqs.  3,  4 and  5,  but  the  value  of 
CpRe  In  Eqs.  1 and  2 had  to  be  calculated  In  each  step  by  numerical 
solution  of  the  definitive  empirical  equations  (Eqs.  6 and  7).  This  was 
done  using  Newton's  Method  (Carnahan  et  al.,  1969)  for  finding  the  zeros 
of  a function.  The  values  of  ux  and  vy  In  Eqs.  1,  2 and  3 were  calculated 
from  Eqs.  8 and  9 for  a cylindrical  target  and  Eqs.  10,  11,  12  and  13  for 
a spherics!  target. 

Once  the  value  of  yd  had  been  determined  the  Integration  pro- 
cedures described  In  Section  8 of  this  report  were  applied.  The  values 
of  f(yj)  used  in  Eqs.  41  and  43  were  obtained  from  Eq.  39.  The  quantities 
v and  y used  In  the  latter  equation  were  calculated  from  Eqs.  23a  and  21a 
respectively.  The  values  of  vx,  vy  and  y'  used  In  the  latter  two  equations 
were  calculated  by  following  a particle  from  Its  Initial  position  to  the 
cylindrical  or  spherical  target  using  the  previously  described  step  by 
step  procedure  applied  In  determining  the  value  of  yd> 

10.  IMPACTION  EFFICIENCY 

The  Impaction  efficiency  was  defined  as  the  ratio  of  the  cross 
sectional  area  of  the  far  upstream  envelope  of  particles  which  eventually 
hit  the  target  to  the  cross  sectional  area  of  the  target  Itself.  For  a 
transverse  cylindrical  target  this  Is  simply  given  by  the  value  of  yd> 

For  a spherical  target  It  Is  given  by  yd2. 

1 1 . COMPUTER  PROGRAMS 

Computer  programs  were  written  In  Fortran  IV  for  the  DRES 
IBM  1130  computer  to  obtain  solutions  to  the  drag  coefficient  equations 
by  the  method  described  In  Section  9 of  this  report.  The  complete  pro- 
grams, together  with  one  set  of  results  for  each,  are  shown  In  Appendix  A 
for  the  cylinder  In  free  space  and  In  Appendix  B for  the  sphere  In  a tube 
or  In  free  space  (Appendix  B).  These  programs  are  annotated  so  that  the 
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functions  of  their  various  parts  can  be  understood  without  further 
description.  The  application  of  the  method  described  In  Section  9, 
which  was  used  In  the  programs,  Is  straight  forward  with  the  exception 
of  the  step  by  step  Integration  procedure  near  the  target  sphere  or 
cylinder.  This  procedure  Is  therefore  explained  as  follows. 

The  particle  motion  was  calculated  In  time  steps,  AT,  until 
the  particle  was  just  one  time  step  away  from  the  target.  Then  the  time 
Increment  size  was  decreased  by  a factor  of  ten  and  step  by  step  Inte- 
gration continued  until  the  particle  was  one  new  time  step  away  from  the 
target.  At  this  point  the  time  Increment  size  was  decreased  by  another 
factor  of  ten  and  the  Integration  allowed  to  proceed  until  the  target 
was  reached.  This  method  ensured  that  the  position  of  the  target  and 
particle  coincided  within  a maximum  error  given  by  the  distance  travelled 
during  one  percent  of  the  original  time  step  size,  while  allowing  the 
particle  to  reach  the  proximity  of  the  target  In  an  adequate  but  small 
number  of  steps. 

12.  RESULTS 

Drag  coefficients  and  impaction  efficiencies  were  calculated 
by  means  of  the  computer  programs  (Appendix  A and  B)  from  two  main  points 
of  view.  First,  they  were  calculated  using  the  particle  sizes,  fluid 
velocities  and  target  configurations  used  In  experiments  done  at  Suffleld 
many  years  ago  (Gillespie  and  Gunter,  1957;  1959).  This  was  done  so  that  a 
comparison  of  the  results  from  the  theoretical  calculation  of  drag  due 
to  particles  In  air  could  be  made  with  existing  Information  obtained 
from  experiments.  The  results  were  obtained  for  a narrow  range  of  particle 
sizes  and  fluid  velocities  and,  correspondingly,  a narrow  range  of  Inertia 
parameters  and  free  stream  particle  Reynolds  numbers.  Next,  calculations 
of  Impaction  efficiencies  and  drag  coefficients  were  done  for  a wide 
range  of  Inertia  parameters.  The  Impaction  efficiencies  were  compared 
with  the  results  of  other  workers  (Frledlander,  1977)  but  no  results  were 
available  for  comparison  to  the  calculated  drag  coefficients  over  this 
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wide  range  of  the  inertia  parameter.  All  the  results  previously  men- 
tioned are  shown  in  Tables  1 and  2 and  Fig.  4 - Fig.  9,  along  with 
plots  of  the  calculated  distribution  of  forces  due  to  elastic  reflection 
of  particles  from  a cylinder  (Fig.  10)  and  a sphere  (Fig.  11).  The 
latter  two  figures  were  plotted  from  the  sample  computer  results  shown 
with  their  corresponding  programs  In  Appendix  A and  B.  Further  details 
explaining  the  significance  of  the  tables  and  figures  are  described  as 
follows. 

The  calculated  results  for  the  various  particle  sizes  and  fluid 
velocities  used  in  the  experiments  done  at  Suffleld  (Gillespie  and 
Gunter,  1957;  1959)  are  shown  for  a cylinder  (Table  1)  and  for  a sphere 
(Table  2).  The  sizes  of  these  two  targets  also  correspond  to  those  In 
the  experiments.  The  calculations  were  done  for  a 0.2  centimetre  diameter 
cylinder  In  free  space  and  a 0.8  centimetre  sphere  In  a 2.54  centimetre 
tube.  The  tube  which  was  used  In  all  the  experiments  was  Included  for 
the  sphere  since  the  velocity  potential  was  readily  available  (Smythe, 
1964).  However,  the  effect  of  the  presence  of  the  tube  on  the  potential 
flow  and  the  final  theoretical  results  was  found  to  be  negligible.  Simi- 
lar 1 ly  , the  effect  of  not  Including  the  tube  In  the  potential  flow  for 
the  cylinder  Is  expected  to  be  negligible., 

A reason  for  tabulating  the  calculated  results  Is  to  show 
clearly  that  the  effect  of  particle  size  on  the  calculated  drag  co- 
efficient for  the  sizes  used  In  the  experiments  was  negligibly  small. 

This  conclusion  was  also  drawn  from  the  axperlments  reported  by  Gillespie 
and  Gunter  (1957;  1959).  The  drag  coefficients  due  to  the  middle  size 
particles  alone  were  plotted  against  target  Reynolds  number  for  the 
cylinder  (Fig.  4)  and  for  the  sphere  (Fig.  5).  Also  shown  for  the  pur- 
pose of  comparison  of  experiment  to  theory,  described  in  the  following 
section,  are  the  drag  coefficients  due  to  air  alone. 

The  Impaction  efficiencies  for  a wide  range  of  the  Inertia 
parameter  were  plotted  for  cylinders  (Fig.  6)  and  for  spheres  (Fig.  8). 

The  cylinder  results  were  plotted  for  $ = 1000.  The  definition  of  4 
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and  the  reason  for  Its  value  choice  are  described  as  follows.  The 
collection  efficiency  and  also  the  drag  coefficient  are  a function  of 
two  dimensionless  groups,  the  Inertial  parameter,  K,  and  the  free  stream 
particle  Reynolds  number,  ReQ.  A new  dimensionless  group 

Re  * 

♦ “ — (Eq.  44) 

K 

Independent  of  particle  size  was  Introduced.  According  to  the  rules  of 
dimensionless  analysis  this  Is  permissible,  but  the  efficiency  Is  still 
determined  by  two  groups  which  for  convenience  were  chosen  to  be  K and 
The  collection  efficiencies  were  calculated  by  means  of  a mechanical 
analog  and  plotted  for  values  of  <j>  by  other  workers  (Brun  et  al.,  1953). 
The  results  for  the  middle  value  of  <f>  from  this  work  are  also  shown  (Fig. 
6)  for  comparison.  The  sphere  results  were  plotted  (Fig.  8)  for  an  Re0 
value  of  128.  Again  this  choice  was  made  for  convenience  of  comparison 
to  other  work.  The  choice  was  the  middle  value  of  seven  for  which  curves 
calculated  from  Invlscld  flow  theory  by  Dorsch,  Saper  and  Kadow  are 
shown  (Frledlander,  1977).  The  curve  Is  also  shown  In  Fig.  8. 

The  calculated  drag  coefficient  due  to  particles  alone  were 
plotted  for  a wide  range  of  Inertia  parameters  for  cylinders  (Fig.  7) 
and  for  spheres  (Fig.  9).  In  these  two  figures  each  curve  was  fitted 
by  eye  through  several  data  points  obtained  by  applying  the  computer  pro- 
grams (Appendix  A and  B)  to  various  particle  diameters  for  the  cylinder 
case  and  various  target  diameters  for  the  sphere  case. 

10.  . COMPARISON  OF  THEORETICAL  RESULTS  TO  EXPERIMENTAL  DATA 

The  experimental  work  previously  done  at  Suffleld  (Gillespie 
and  Gunter,  1957;  1959)  Indicated  that  the  drag  force  of  a sphere  or 
cylinder  Is  given  by  the  following  equation: 

F - l/2p*U2CD  A (Eq.  45) 

e 

where  p*  * p + pp  (Eq.  46) 
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and  generally  p and  pp  are  unequal.  For  Eq.  45  to  be  true,  the  drag 
coefficients  due  to  particles  alone  carried  by  air  and  due  to  air  alone 
must  be  equal.  This  Is  shown  as  follows.  The  total  drag  force  con- 
sisted of  two  parts  expressed  by  the  following  equation: 

K + K - 1/2pnU2ACn„  + l/2pU2ACn  (Eq.  47) 


Fp  + Fa  - l/2ppU2ACDp  + l/2pU2ACD  (Eq.  47) 

Comparison  of  Eq.  45  and  46  to  Eq.  47  shows  that 

(p  ♦ Pp)CDe  ■ PpCDp  + PCD  (Eq.  48) 

CDp  and  CD  must  be  equal  for  Eq.  48  to  be  true. 

The  results  obtained  herein  (Tables  1 and  2)  and  Illustrated 
In  Fig.  4 and  5 show  that  the  calculated  drag  coefficients  due  to  particles 
carried  by  air  are  considerably  greater  than  the  established  drag  coef- 
ficients due  to  air  alone  (Schllctlng,  1960).  The  discrepancies  between 
the  experimental  and  calculated  drag  coefficients  were  accounted  for  by 
showing  that  the  velocity  of  the  particles  In  the  experiments  were  con- 
siderably lower  than  the  mean  air  velocity. 

The  particles  which  were  Introduced  by  Gillespie  and  Gunter 
Into  the  air  stream  with  zero  velocity  In  the  direction  of  air  flow  were 
accelerated  along  a horizontal  section  of  2.54  cm  tube,  17  cm  long,  and 
a vertical  section,  78  cm  long,  giving  a total  tube  length  of  95  cm  between 
the  starting  point  and  the  test  section.  Calculations  of  velocity  as  a 
function  of  distance  for  particles  accelerated  from  rest  were  done  to 
simulate  Gillespie's  experiments  (Mellsen,  In  draft).  The  results  for  a 
travelled  distance  of  95  cm  for  the  various  mean  air  velocities  and  number 
median  diameter  used  In  the  experiments  are  shown  In  Table  3.  Also,  the 
curves  of  velocity  versus  distance  for  the  three  diameters  and  one  air 
velocity  are  Illustrated  In  Fig.  14. 

The  particle  velocities  were  difficult  to  calculate  because 
the  behavior  of  the  particles  travelling  In  the  air  stream  around  the 
corner  In  the  tube  was  very  complex.  The  velocity  of  the  air  In  a 
straight  cross  section  of  the  tube  was  zero  at  the  walls  and  Increased 
to  the  centre  to  1.22  to  1.25  times  the  mean  velocity  (Prandtl  and  Tletjens, 
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1957).  Therefore,  the  particle  velocities  were  affected  by  the  location 
of  the  particles  In  the  tube  cross  section.  This  could  not  be  accounted 
for  because  the  time  history  of  the  location  of  the  particles  was  not 
known.  As  a simplifying  approximation,  the  mean  air  velocity  was  assumed 
to  act  on  the  particles  over  their  length  of  travel  In  the  tube  wh^ch  was 
also  assumed  vertical  over  Its  entire  length.  The  experimental  results 
Indicate  that  the  variation  In  velocities  was  less  than  this  method  of 
calculation  shows.  A mean  value  of  drag  force  due  to  the  presence  of 
all  particle  sizes  was  used  to  offset  this. 

The  drag  coefficients  which  were  calculated  for  particles  with 
a free  stream  velocity  equal  to  the  mean  air  velocity  (Table  1 and  2) 
were  then  adjusted  to  simulate  the  special  conditions  of  the  experiments 
as  follows.  The  drag  coefficients  In  the  experiments  were  based  on  the 
mean  air  velocities  as  measured  by  a rotameter.  The  drag  coefficients  due 
to  the  particles  travelling  at  lower  velocities  were  adjusted  to  the  mean 
air  velocity  by  the  following  relationship  for  flow  over  the  cylinder. 


v \2 
rms  1 


(Eq.  49) 


This  method  was  justified  because,  for  the  purpose  of  comparison,  the 
calculated  drag  coefficients  over  the  range  of  velocities  and  particle 
sizes  used  were  sufficiently  near  a constant  value  (Table  1).  All  three 
particle  sizes  were  assumed  present  in  equal  quantities  In  the  air  flow. 
Then,  since  the  drag  forces  vary  as  the  square  of  the  velocity,  the 
root  mean  square  particle  velocity,  vrms,  was  calculated  for  each  mean 
fluid  velocity  (Table  3)  and  used  In  Eq.  49  to  obtain  a value  for  the 
drag  coefficient  due  to  dust  alone  at  each  mean  air  velocity  (Table  4). 
The  overall  drag  coefficients  due  to  air  and  particles  were  then  cal- 
culated for  the  dust  concentrations  used  In  the  experiments,  namely  2 
and  5 kilograms  per  cubic  metre  of  space.  This  was  done  by  solving 
Eq.  48  for  CQ  . That  Is 


CD,  " 


pdC0c  + »CD 


P + P, 


(Eq.  50) 
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The  results  are  shown  for  each  mean  air  velocity  (Table  4)  and  plotted 
along  with  the  drag  coefficient  for  air  alone  (Fig.  12)  as  was  done  for 
the  experiments  of  Gillespie  and  Gunter  (1957;  1959). 


A similar  method  was  used  to  obtain  the  overall  drag  coef- 
ficients for  the  sphere  with  the  exception  that,  as  for  the  experimental 
results,  they  were  corrected  for  the  presence  of  the  tube  walls  by  the 
Landenburg  correction  (Gillespie  and  Gunter,  1957).  That  Is,  the  drag 
coefficients  due  to  particles  alone  were  divided  by  the  factor  1 + 2.4 
L/a  to  obtain  the  corrected  drag  coefficient,  C^.  For  the  sphere  and 
tube  radii  used,  the  factor  Is  1.756.  Then,  the  overall  drag  coefficient 
for  the  sphere  was  obtained  from 


'De 


PPCDL  * pCC 
P + P* 


(Eq.  51) 


The  results  are  shown  for  each  mean  air  velocity  in  Table  4,  and  plotted 
along  with  the  drag  coefficient  for  air  alone  In  Fig.  13. 


14.  DISCUSSION 

The  equations  of  motion  (Eq.  1 and  2)  have  never  been  proven 
experimentally  (Frledlander,  1977)  although  they  have  been  used  to  cal- 
culate deposition  efficiencies  for  cylinders  (Glauert,  1940).  The 
mathematical  model  described  herein  used  these  equations  to  calculate 
collection  efficiencies  for  cylinders  which  compare  favorably  (Fig.  6) 
with  the  results  obtained  using  a differential  analyzer  (Brun  et  al., 
1953).  Also,  the  Impaction  efficiencies  for  spheres  (Fig.  8)  agree  with 
solutions  to  the  equations  obtained  by  other  workers  (Frledlander,  1977). 

In  the  work  reported  herein,  the  mathematical  model  for  cal- 
culating Impaction  efficiency' was  extended  to  obtain  the  Impaction  forces 
due  to  particles  In  air  flowing  over  spheres  and  cylinders.  These  results 
agree  favorably  with  experiments  previously  done  at  Suffleld  (Gillespie 
and  Gunter,  1957)  when  adjusted  to  the  special  conditions  of  the  experi- 
ment. The  experiments  do  support  the  theory,  but  the  theoretical  results 
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show  that  the  results  of  the  experiments  are  not  applicable  to  the  free 
field  case  of  drag  due  to  dusty  air  flowing  over  a cylinder  or  sphere. 

The  model  could  also  be  extended  to  blast  loading  problems. 

Since  the  effect  of  compressibility  becomes  Important  In  these  problems 
It  would  have  to  be  accounted  for.  This  could  be  done  In  the  mathe- 
matical model  by  changing  the  air  flow  field  equations  to  Include  It. 

The  model  could  also  be  used  to  calculate  the  actual  distributions  of 
Impaction  velocities  on  cylinders  (Fig.  10)  and  spheres  (Fig.  11)  in 
special  cases  which  may  have  applications  In  chemical  and  nuclear  defence. 

15.  CONCLUSIONS 

A mathematical  model  was  developed  for  calculating  the  Impaction 
forces  due  to  airborne  particles  on  spheres  and  cylinders.  The  model  can 
also  be  used  for  calculating  Impaction  efflclences,  velocity  distributions, 
and  points  of  Impact  for  a stream  of  airborne  particles  flowing  over  a 
sphere  or  cylinder.  Hence,  the  drag  coefficient  for  the  particles  on  the 
sphere  or  cylinder  can  be  determined. 
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TABLE  3 

VELOCITY  OF  PARTICLES  AT  TARGET  POSITION 


ACCELERATED  VERTICALLY  UPWARD  IN  A TUBE  FROM  REST 


95  cm  UPSTREAM 


PARTICLE 

MEAN  AIR 

PARTICLE  VELOCITY 

RMS  PARTICLE 

DIAMETER 

VELOCITY  IN 

AT  TARGET 

VELOCITY  AT  TARGET 

cm 

TUBE  cm  sec*1 

POSITION  cm  sec-1 

POSITION  cm  sec"1 

777 


l 299 

1 572 

534 

1 664* 

* Upward  terminal  velocity  reached 
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NOTION  OF  ThF  PARTICLES  at  the  cylinder  position  is  given  by 
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the  «»OTION  OF  A CRITICAL  PARTICLE  IS  GIVEN  BY 
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MOTION  OF  THE  PARTICLES  AT  THE  SPHERE  POSITION  IS  GIVEN  BY 
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